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f-y  ^ 

where  kj^)  = k&aJ  + ^2h)*  measurements  were  made  in  a flow  tube  apparatus 

at  T = 302  {£)2*K  at  a total  pressure  of  3 Torr.  0(^P)  was  produced  by  the 
titration  of  active  nitrogen  with  NO.  H2(v  = 1)  was  produced  by  a heated  tung- 
sten filament,  and  its  concentration  was  d^ermined  from  VUV  absorption  mea- 
surements in  the  Lyman  band  system  of  Absolute  concentrations  of 

OH(v  = 0)  and  OH(v  = 1)  were  deduced  from  laser  induced  fluorescence  in  OH 
in  a calibrated  fluorescence  system,  r— The  rates  and  uncertainties  are; 
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= (‘-otre) 


sec 
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X 10  cc /molecule -sec 
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^ 4.7X10"  cc/molecule -sec 


The  reported  value  for  represents  an  increase  owing  to  vibrational  excitation 
in  hydrogen  by  a factor  of  <=»2.6  X 10^  over  the  fully  ground  state  reaction  rate 
constant.  Measurements  of  the  fully  ground  state  reaction  rate  constant  and 
recent  theoretical  treatments  of  the  O -t  H2  reaction  are  discussed  in  relation 
to  these  results. 
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I.  INTRODUCTION 


The  reaction  of  0(^P)  with  has  been  the  subject  of  many  investi- 

1-4 

gations  in  the  last  twenty  years.  A review  of  all  of  the  measurements 
through  June  1971  is  presented  in  Ref.  5.  In  the  temperature  range  from 

O O 

400  K to  2000  K there  is  general  agreement  as  to  the  value  of  the  rate 
constant.  There  have  been  no  reported  measurements  of  the  rate  of  this 
reaction  for  vibrationally  excited  hydrogen.  The  fully  ground  state  reac- 


0(^P)  + = 0)— ►OH(v  = 0)  + H 

is  endothermic  by  about  2.  0 kcal/mole  and  has  an  activation  energy  of 

5 

about  8.  9 kcal/mole.  When  the  H-,  has  one  quantum  of  vibrational 
energy,  reaction  (2a) 

0(^P)  + = 1)— ►OH{v  = 1)  + H 

is  barely  endothermic  («  0.3  kcal/mole)  and  reaction  (2b) 

0(^P)  + H2(v  = 1)— ►OH(v  = 0)  + H 

is  exothermic  by  9.  9 kcal/mole.  Thus,  it  is  seen  that  if  vibrational 

energy  in  is  effective  in  overcoming  the  activation  energy  of  reaction 

(1),  the  ratio  k^/k^,  where  k2  = k2^  + k2^,  might  be  extremely  large. 

Birely  et  al.  ^ have  reported  an  upper  limit  to  the  value  of  k2  such  that 
4 ^ 

k_/k,  s 3,4  X 10  at  T = 300  K.  This  means  that  at  most  the  fraction 

a = 0.  52  of  the  H2  vibrational  energy  is  effective  in  reducing  the  activa- 
■7 

tion  energy.  In  addition,  there  have  been  two  reports  of  theoretical 
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treatments  of  this  reaction.  ' Johnson  and  Winter  report  an  Arrhenius 


expression  for 

k2  = 4.  65T  X lO"^^  exp(-1868/T)  cc/molecule-sec. , 

which  is  lower  than  the  upper  bound  established  by  Birely  et  al.  at 
T = 300  “k. 

Interest  in  the  dependence  of  elementary  reaction  rate  constants  on 
particular  quantum  states  of  reactants  and  products  has  grown  enormously 
in  the  past  two  decades.  These  data  are  important  in  characterizing  the 
detailed  mechanism  and  potential  energy  hypersurface  for  a reaction,  and 
in  more  practical  application  to  such  areas  as  molecular  lasers,  laser- 
initiated  chem.istry,  laser  isotope  separation  and  chemical  processes 
occurring  under  a variety  of  other  nonequilibrium  situations.  Molecular 
beam,  chemiluminescence  and  chemical  laser  studies  have  led  to  the 
characterization  of  energy  partitioning  in  a variety  of  elementary  reac- 
tions. However,  the  effect  of  reactant  vibrational  energy  is  not  well 
characterized  and  measurements  which  are  differential  with  respect  to 
both  reactant  and  product  quantum  states  are  very  rare. 


y 


! 
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This  paper  reports  for  the  first  time  a measurement  of  the  rate  of 
the  reaction 


3 

OC’P)  + H2(v  = 1)  — ► OH(v  = 1)  + H. 
Also  reported  is  an  upper  limit  to  the  value  of  where 


(2a) 


II.  EXPERIIVLENTAL 
a.  Apparatus 

The  measurements  were  conducted  at  302  ± 2 K in  the  flow  tube 

3 

apparatus  shown  schematically  in  Fig.  1.  Atomic  oxygen,  0(  P),  was 
generated  by  the  rapid  reaction  N + NO — ► + O,  for  which  the  N-atoms 

were  provided  by  a microwave  discharge  in  pure  N^.  Nitric  oxide  was 
admitted  to  the  flow  tube  through  a circular  ring  of  small  radial  holes  In 
the  wall  of  the  flow  tube.  Bebween  the  N^  microwave  discharge  and  the 
NO  titration  inlet  there  was  placed  a side  arm  for  the  admission  of  bypass 
N2  gas.  The  purpose  for  this  bypass  flow  is  discussed  later.  Some  14  cm 
downstream  from  the  NO  titration  inlet  the  O-atoms  were  allowed  to  mix 
with  pure  H gas  containing  vibrationally  excited  H.  (v  =1).  Vibrational 
excitation  in  the  hydrogen  was  created  by  passing  the  through  a 1 cm 
diameter  alumina  tube  coi.  taining  a densely  wound  tungsten  filament  which 
was  electrically  heated.  Upon  exiting  from  this  gas  heater,  rotational 
and  translational  modes  of  the  were  equilibrated  with  the  room  temper- 
ature flow  tube  walls  before  mixing  with  the  O-atom  stream.  However, 
the  VUV  absorption  measurements  confirmed  that  tlie  vibrational  modes  of 
H_  equilibrate  much  more  slowly  with  the  walls.  Hence,  we  have  a slowly 
decaying  H_(v  = 1)  density  as  a function  of  distance  traveled  down  the  flow 
tube.  Auxiliary  measurements  without  the  O-atom  source  confirmed  the 
absence  of  any  OH  production  by  a possible  oxidation  of  the  H^  on  the  hot 
alumina  surfaces.  An  annular  cooling  water  jacket  was  provided  in  the 
flow  tube  walls  between  the  bottom  of  the  gas  heater  and  the  reactant  mix- 
ing point  for  the  purpose  of  ensuring  a controlled  wall  temperature. 


IL 
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From  the  point  where  the  O and  reactant  streams  meet  they 
could  be  directed  into  either  of  two  portions  of  the  flow  tube.  They  could 
be  convected  into  the  VUV  absorption  cell  where  quantitative  measure- 
ments of  the  density  of  = 0)  and  = 1)  were  performed.  Con- 

versely, the  gas  flow  could  be  convected  to  a section  of  flow  tube  designed 
for  the  generation  and  collection  of  laser  induced  fluorescence  in  the  pro- 
duct OH  molecules.  The  OH  fluorescence  cell  and  the  center  of  the  VUV 
absorption  cell  were  equidistant  from  the  point  where  the  reactant 
streams  mixed.  The  light  source  for  the  VUV  absorption  measurements  was 
a pulsed  dc  electrical  discharge  in  argon.  This  discharge  produces  the  argon 
continuum  between  110  cuid  140  nm.  This  emission  passes  through  a lithium 
fluoride  (UiF)  window  into  a 14-cm-long  absorption  path  in  the  flow  tube.  It 
then  exits  from  the  flow  tube,  passes  through  yet  another  LiF  window  eind 
through  the  entrance  slit  of  a sceuining  VUV  monochromator.  In  operation,  the 
monochromator  is  scanned  slowly  through  a particular  rotation -vibration  line 
in  the  Lyman  band  system  of  H^*  For  measurements  of  H^  (v  = 1)  density, 
absorption  by  the  P(l),  R(2)  line  blend  of  the  (1,  1)  band  at  \ = 114.62  nm  was 
measured.  The  resulting  absorption  trace  was  manually  integrated  to  yield 
equivalent  width.  The  density  of  H2(v  = 1)  was  then  deduced  from  curve -of- 
growth  calculations.  The  accuracy  of  the  curve -of-growth  analysis  was 
verified  by  absorption  measurements  in  the  (0,  0)  band  of  this  system  and  cor- 
responding pressure  measurements  with  a bourdon  type  pressure  gauge. 

For  the  laser  induced  fluorescence  measurements  of  OH  density,  focused, 
tunable  laser  radiation  entered  the  flow  tube  through  a quartz  window  oriented 
at  Brewster's  angle.  It  then  passed  through  some  optical  baffles,  was  directed 
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along  the  flow  tube  axis  through  more  optical  bciftles  and  exited  through  another 
quartz  Brewster  angle  window.  Laser  pulse  energy  was  measured  by  a pyro- 

V electric  detector  in  the  exit  beam.  Laser  induced  fluorescence  was  collected 

^ by  a quartz  lens  on  one  side  of  the  flow  tube  and  by  an  aluminized  spherical 

mirror  diametrically  opposite  the  lens.  These  two  elements  were  arreinged 
, > with  an  optical  aperture  of  f/1. 5.  Collected  radiation  then  passed  through  a 

y broadband  interference  filter  and  directly  into  a photomultiplier  tube 

V 

(EMI  9635QB).  The  filter  had  a peak  transmission  of  Z4%  at  \ = 312.  5 nm, 

12%  at  \ = 305.  5 nm  and  X.  = 318.0,  and  less  than  1%  for  \ < 299.5  nm  and 
X.  > 328.  5 nm.  Hence,  this  filter  readily  passed  radiation  in  the  (0-0)  and 
“ * (1-1)  bands  of  the  OH(A  2 —►X  IT)  system.  The  PMT  anode  was  terminated 

in  a 50  ohm  resistor  at  the  input  of  a preamplifier  (70  MHz  bandpass,  gain  = 10). 
The  preamp  output  was  connected  to  a PAR  Model  162  Dual  Ch«innel  Boxcar 
Averager. 


The  excLtation  source  was  a Molectron  DL-200  tunable  dye  laser. 

The  laser  output  was  frequency  doubled  with  a 4 cm  long  angle  tuned 

crystal  of  ADP.  After  frequency  doubling  to  X = 308.  5 nmi  the  typical 

laser  pulse  contained  5 to  10  X 10  ^ joules  of  energy,  was  approximately 
-9  -1 

5 X 10  sec  in  duration,  and  1.0  cm  in  spectral  width.  The  laser  was 
operated  at  a pulse  rate  of  25  Hz. 

The  flow  tube  was  constructed  of  Pyrex  glass,  with  an  internal  dia- 
meter of  25  n.m.  The  average  gas  flow  velocity  was  15  m/sec.  Exper- 
iments were  conducted  with  a fixed  total  pressure  of  3.0  Torr  comprised 
of  equal  parts  of  and  H^.  The  exterior  surface  of  the  flow  tube  was 


reactants  mixed,  one  at  the  fluorescence  cell  and  the  third  midway 
between  these  two.  Most  of  the  Inside  surface  of  the  flow  tube  was  coated 
with  a halocarbon  wax  which  has  been  shown^^  to  be  effective  in  control- 
ling the  heterogeneous  removal  rate  of  various  atoms  and  radicals  at  the 
walls.  No  coating  was  ^.pplied  in  the  region  near  the  hot  tungsten  filament. 
In  the  immediate  vicinity  of  the  microwave  discharge  boric  acid  was 
applied  to  the  tube  walls.  The  acid  was  converted  to  the  powdery  white 
oxide  by  heating  the  glass  from  the  outside  with  a propane  torch. 

The  flow  rates  of  and  gases  were  monitored  by  Matheson 
rotameter  gauges.  The  flow  rate  of  nitric  oxide  was  determined  by  meas- 
uring the  pressure  drop  in  a known  reservoir  volume.  The  O-atom  den- 

12 

sity  was  deduced  at  the  visible  chemiluminescence  titration  end  point 
from  the  known  NO  flow  rate  and  computer  integration  of  a set  of  twelve 
chemical  reactions.  These  same  calculations  provided  the  densities  of 
several  minor  species  which  accompany  the  N + NO  system  and  whose 
role  in  the  overall  kinetics  of  OH  production  and  removal  must  be  under- 
stood. 

The  gases  used  in  the  experiment  all  were  supplied  by  Matheson 
Company  with  the  following  stated  purities:  UHP  grade  99.  999%  min.  , 

N^,  UHP  grade  99.  999%  min. , and  NO,  C.  P.  grade,  99.0%  min.  Both 
the  H^  and  N^  were  further  purified  by  passing  through  traps  filled  with 
glass  beads  at  77 "k  and  700  Torr  pressure.  The  nitric  oxide  was  puri- 

O O 

fled  by  three  successive  distillations  between  traps  at  77  K and  156  K 
(ethanol  slush). 
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b.  Procedure  and  Calibration 

Fluorescence  amplitude  data  were  obtained  by  tuning  the  laser  out- 
put  frequency  through  a particular  rotational  line  in  the  (A  E — ►X  H.) 
system  of  OH.  This  procedure  provided  a measure  of  the  background  sig- 
nal and  the  sum  of  the  background  and  fluorescence  signals  for  each 
observation.  The  Qj^(l)  line  was  used  in  both  the  (0-0)  and  the  (1-1)  bands 
to  collect  detailed  data  on  OH(v)  density.  Fluorescence  radiation  was  col- 
lected from  the  (0-0)  band  when  the  v"  = 0 state  was  investigated  and  from 
both  the  (0-0)  and  the  (1-1)  bands  simultaneously,  when  the  v"  = 1 state 
was  investigated. 

The  density  of  nitric  oxide  and  of  several  other  important  species 
was  conveniently  varied  by  changing  the  ratio  of  nitric  oxide  inlet  flow 
rate  to  N-atom  flow  rate.  This  ratio  is  hereafter  referred  to  as, 

I = ® stands  for  flow  rate  in  units  of  particles/sec.  The 

computed  dependences  of  the  densities  of  O,  N,  NO,  and  NO^  on  § are 
shown  for  a typical  case  in  Fig.  2.  The  kinetic  scheme  used  to  compute 
these  results  is  detailed  in  Appendix  A.  As  a test  of  the  ability  of  the 
computer  code  for  numerical  integration  of  the  chemical  rate  equations 
to  accurately  simulate  the  chemistry  in  this  experiment  it  was  decided  to 
measure  the  density  of  NO2  produced  by  the  N + NO  system  for  compari- 
son with  calculated  values.  For  this  purpose  calibrated  laser  induced 
fluorescence  data  were  obtained  for  NO^.  The  laser  was  tuned  to  445.  0 nm 
and  fluorescence  was  observed  with  a photomultiplier /filter  comMnation 
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Fig.  Z.  Calculated  Species  Concentrations  As  a Fvinction  of 
NO  Flow  Rate  for  the  Reactions  Listed  in 
Appendix  A 
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for  X s 480  nm.  Data  were  obtained  for  fluorescence  Intensity  as  a func- 
tion of  § for  one  value  of  total  N-atom  flux.  The  data  and  corresponding 
computer  predictions  are  shown  in  Fig.  3.  It  can  be  seen  that  the  agree- 
ment is  good.  One  therefore  has  a measure  of  confidence  in  the  charac- 
terization of  the  kinetics  of  this  system  and  in  the  validity  of  the  computer 
calculations. 

The  density  of  ^^^(v  = 1)  could  be  readily  varied  by  changing  the 
temperature  of  the  tungsten  filament.  However,  it  was  discovered  early 
in  the  investigation  that  vibrationally  excited  hydrogen  was  sometimes 
present  in  the  flow  even  when  the  tungsten  filament  was  turned  off,  if  the 
microwave  discharge  in  N2  was  on.  This  effect  was  unimportant  for  any 
of  the  data  presented  in  this  paper.  The  exact  nature  of  the  energy  trans- 
fer process  is  not  known  but  the  available  evidence  suggests  a V-V  mech- 
anism between  vibrationally  excited  N^,  with  at  least  two  quanta  of  vibra- 
tional energy,  and  the  cold  hydrogen  gas.  Very  high  levels  of  vibrational 
excitation  are  known  to  be  produced  by  a microwave  discharge  in  N^. 

It  was  further  found  that  there  were  two  methods  by  which  this 
energy  transfer  was  prevented.  By  adding  a small  amount  of  CO^  gas  to 
the  discharged  N2  stream  before  it  contacted  the  flow,  the  generation 
of  i 1)  by  this  process  was  prevented.  It  was  also  verified  that  there 

was  no  effect  on  the  density  of  H2(v  = 1)  if  CO^  was  titrated  into  the  flow 
in  the  absence  of  active  nitrogen.  The  CO2  evidently  quenches  the  active 
nitrogen  witliout  directly  affecting  H2{v  = 1).  The  second  method  involved 
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adding  an  excess  of  undischarged  downstream  of  the  microwave  dis- 
charge; the  generation  of  vibrationally  excited  by  this  unknown  energy 
transfer  process  was  prevented  without  any  effect  on  H2(v  =1)  coming 
from  the  tungsten  filament. 

Although  this  unknown  energy  transfer  process  generates  copious 
amounts  of  H2(v  =1),  more  than  can  be  produced  by  the  tungsten  filament 
alone,  it  was  deemed  unadvisable  to  use  it  as  a source  of  reactant  in  the 
present  experiments.  All  of  the  data  presented  in  this  paper,  therefore, 
were  obtained  when  the  amount  of  = 1)  produced  by  this  unknown  pro- 

cess was  undetectable.  For  most  of  the  data  presented,  this  meant  that 
the  H2(v  = 1)  density  exceeded  that  due  to  this  unknown  energy  transfer 
process  by  at  least  a factor  of  30.  The  bypass  flow  method  of  prevent- 
ing this  energy  transfer  process  was  used  exclusively  in  these  data.  The 
total  flow  rate  was  held  constant  while  the  relative  proportion  of  bypass 
flow  was  varied. 

The  possible  removal  of  reactive  species  by  heterogeneous  reactions 

at  the  wall  must  be  understood  in  order  to  properly  interpret  the  results. 

As  noted  in  Section  (a)  above,  the  walls  throughout  the  reaction  zone  are 

coated  with  a fully  halogenated  fluorocarbon  wax.  As  reported  originally 

by  Glass,  this  coating  is  effective  in  preventing  H-atom  and  O-atom 

13 

recombination  on  the  walls.  Glass  also  has  reported  a wall  removal 
rate  expressed  as  a first  order  rate  constant  = -30  sec  ^ for  OH  and 

a vibrational  relaxation  rate  for  OH(v  = 1)  ^ -80  sec~^  when  such  a 
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wax  Is  used  to  coat  the  walls.  These  rates  are  sufficiently  slow  that  the 
removal  of  these  species  is  dominated  by  homogeneous  gas  phase  reac> 

tions  under  all  conditions  encountered  in  ttiis  experiment.  For  H (v  = 1), 

6 ^ 

Bwrely  et  al.  reported  in  the  range  -62  to  -91  sec"^  for  an  uncoated 

Pyrex  wall.  We  have  measured  ap  average  value  of  = -46  sec'^  for 
H^Cv  = 1)  deactivation  on  a halocarbon  wax  coated  wall.  This  value  was 
derived  from  measurements  of  the  density  of  H^iv  = 1)  as  a function  of 
the  exposed  length  of  coated  glass  tube.  Because  H2(v  = 1)  does  not  react 
rapidly  with  any  species  present  in  these  experiments,  this  wall  deactiva- 
tion has  a minor  effect  on  the  OH(v)  kinetics.  As  will  be  documented 
later,  the  homogeneous  gas  phase  removal  processes  for  both  OH(0)  and 
OH(l)  are  sc  fast  that  it  is  of  primary  importance  to  know  the  density  of 
H2(v  = 1)  only  at  the  point  where  measurements  of  OH(v)  densities  are 
made  and  of  only  secondary  importance  to  know  the  time  history  of 
H2(v  = 1)  density  along  the  tube. 

Wavelength  calibration  of  the  exciting  laser  radiation  was  accom- 
plished by  obtaining  an  excitation  spectrum  for  OH  laser  induced  fluores- 
cence on  a propane  flame  at  atmospheric  pressure.  Fluorescence  was 
detected  from  more  than  one  hundred  separate  lines  in  the  (0-0)  and  (1-1) 
bands  of  OH  with  a standard  deviation  of  0.  07  A between  measured  wave- 
length and  wavelengths  for  OH  taken  from  Dieke  and  Crosswhite. 

As  will  be  documented  later,  the  loss  processes  for  both  OH(v  = 0) 
and  OH(v  = 1)  are  sufficiently  rapid  and  the  depletion  of  reactants  is 


r- 

> 
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sufficiently  slow  that  a steady  state  analysis  of  the  data  is  required  where- 
in it  is  necessary  to  obtain  the  absolute  density  of  OH(v)  in  order  to  obtain 
rate  constants. 

The  reaction  H + NO^ — ►OH(v)  + NO  was  used  as  a quantifiable 

source  of  OH(v  = 0)  in  order  to  calibrate  the  fluorescence  amplitude. 

Details  of  this  procedure  with  corroborating  experimental  evidence  aind  a 

15 

discussion  of  the  underlying  chemical  kinetics  will  be  given  elsewhere. 

Briefly,  the  H + NO^  reaction  is  relied  on  only  as  a source  of  OH(v  = 0). 

The  calibration  for  OH(v  = 1)  comes  from  supplementary  measurements 

in  the  O + system.  The  justification  for  using  H + NO^  rests  partly  on 

13 

the  recently  published  values  for  the  rate  of  physiccd  deactivation  of 
OH(v)  by  H and  by  NO.  These  rates  are  sufficiently  rapid  that  the  resul- 
tant density  of  OH(v  = 0)  after  about  10  msec  of  reaction  time,  is  reason- 
ably well  defined  irrespective  of  the  present  uncertainty  in  the  branching 

13  1 6 17 

ratios  for  producing  different  OH  vibrational  states.  ' ' It  is  true 

that  the  magnitude  of  the  quenching  rates  of  OH(v)  by  H and  by  NO  as 
13 

reported  by  Glass  are  dependent  on  a knowledge  of  just  these  branching 
ratios.  However,  the  present  uncertainty  in  these  branching  ratio  values 
and  the  corresponding  uncertainties  in  the  vibrational  quenching  rates  do 
not  appreciably  influence  the  calculated  OH(v  = 0)  densities. 


1 


The  absolute  sensitivity  for  detection  of  OH(v  = 0)  in  the  H + NO^ 
system  was  obtained  for  laser  excitation  on  the  Q^(l)  line  of  the  (0,  0)  band. 
In  order  to  determine  the  calibration  in  the  O + H^  system  a correction 
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2 + 

must  be  made  for  differences  in  electronic  quenching  of  OH(A  Z , v'  = 0) 
between  the  two  gas  systems.  The  only  significant  difference  is  in  the 
pressure  of  which  is  1.  5 Torr  in  the  case  of  O + H2  and  is  appreciably 
lower  in  the  case  of  the  H + NO2  data.  Measurements  were  obtained  for 
the  quenching  coefficient  in  N^,  and  in  and  an  upper  limit  was  obtained 
in  argon  diluent.  The  fluorescence  decay  rate  was  determined  by  slowly 
scanning  a narrow  gate  (12  nsec  aperture  duration)  across  the  fluores- 
cence signal  using  the  Boxcar  Averager.  An  example  of  the  fluorescence 
decay  data  obtained  is  shown  in  Fig.  4.  Data  of  this  type  were  obtained  at 
a number  of  different  values  of  quenching  gas  pressure.  The  deduced 

quenching  coefficients  are  shown  in  Table  I along  with  some  recently  pub- 
18 

lished  values.  The  agreement  is  seen  to  be  remarkably  good. 

The  relevant  expression  for  the  density  of  OH(v  = 0)  is: 

[OH(v  = 0)]  = K (■^)  . ( 

\'^rs/(o-o) 

where 

Kq  is  the  calibration  coefficient  for  OH(v  = 0),  valid  for  this  system 
only. 

V is  the  peak  signal  voltage  as  the  laser  is  scanned  through  the 
8 

Qj(l)  line. 

Q|  is  the  laser  pulse  energy. 

f is  a dimensionless  factor  which  corrects  for  nonlinear  optical 
s 

absorption,  i.  e. , partial  optical  saturation. 
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RELATIVE  FLUORESCENCE  AMPLITUDE 


2 + 

Table  I.  Electronic  Quenching  Rates  of  OH(A  E , v'  = 0) 
(in  cc/molecule-sec) 


Quenchant 

k 

q 

k 

q 

(Ref.  18) 

(This  work) 

»2 

1.3  X 10"^° 

1.  35(±  . 14)  X 10"^® 

^2 

1. 95  X 10"^^ 

1.  99(±  .49)  X 10"^^ 

Ar 

-13 

5 X 10  ^ 

< 1. 5 X 10"^^ 
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was  observed  in  the  OH  fluorescence 


1 9-22 

Nonlinear  optical  absorption 
data  collected  in  this  experiment.  Data  which  clearly  show  the  transition 
from  linear  to  nonlinear  behavior  are  shown  in  Fig.  5.  These  were 
obtained  by  placing  filters  of  successively  greater  optical  density  in  the 
laser  beam  before  it  entered  the  fluorescence  cell.  The  linearity  of  the 
photomultiplier  detection  system  was  established  in  separate  experiments. 

The  functional  form  to  be  expected  for  the  dependence  of  fluorescence 
on  incident  laser  beam  energy  would  be  difficult  to  calculate  here.  There 
is  a wide  distribution  of  power  densities  represented  in  the  field  of  view  of 
the  photomultiplier  detector  because  the  laser  beam  is  focused  in  the  cen- 
ter of  this  field.  The  observed  nonlinearity  is  effectively  the  result  of  the 
volume  integral  of  the  fluorescence  over  the  irradiated  region  as  modified 
by  the  collection  efficiency  of  the  optics  for  each  point  in  the  field.  In  addi- 
tion, it  may  be  necessary  to  account  for  power  broadening  of  the  absorp- 

23 

tion  line  as  done  by  Killinger  et  al.  Our  approach  is  strictly  empirical. 

The  factor  f in  Eq.  (a)  is  the  ratio  of  the  observed  fluorescence  signal  to 
s 

what  the  signal  would  have  been  if  the  nonlinearity  had  not  set  in.  It  was 
determined  from  data  such  as  are  shown  in  Fig.  5.  Similar  data  were 
obtained  for  every  transition  used  for  quantitative  fluorescence  measure- 
ments. 

A somewhat  different  procedure  had  to  be  used  to  obtain  a calibra- 
tion for  OH(v  = 1)  since  quantifiable  sources  are  not  readily  available. 

With  a known  density  of  OH(X^n^,  v"  = 0,  N"  = 1,  J"  = 3/2),  i.  e.,  the 


I 


RELATIVE 

FLUORESCENCE 

AMPLITUDE 


RELATIVE  LASER  ENERGY 


Experimental  Data  That  Show  Nonlinear  Fluorescence  Behavior, 
Owing  to  Partial  Optical  Saturation.  The  laser  excitation  was 
on  the  Q^(l)  line  of  the  (0,  0)  band  of  OH.  The  laser  beam  energy 
was  5.4  pj  at  a value  of  1.0  on  the  relative  energy  scale. 


lower  level  of  the  Qj(l)  transition,  the  laser  was  tuned  to  the  0^(1)  line 
in  the  (v'  = 1,  v"  = 0)  band  at  281.  915  nm.  The  lower  level  of  this  transi- 
tion is  the  same  as  that  for  the  Qj(l)  line  in  the  (0,0)  band.  The  upper 
level  of  this  transition  is  the  same  as  the  upper  level  in  a Qj(l)  transition 
in  the  (1,  1)  band.  The  fluorescence  intensity  observed  for  (1,0)  band 
excitation  differs  from  what  would  have  been  observed  for  (1,  1)  band  exci- 
tation only  by  the  differing  densities  of  the  lower  states,  i.  e.  , the  ratio 
[OH(v"  = l)]/[OH(v"  = 0)],  and  by  the  known  transition  probabilities.  A,  and 
wavelengths,  laser  pulse  energies,  and  nonlinear  absorption  correction 
factors.  No  corrections  are  necessary  for  differences  in  fluorescence 
spectrum  because  the  upper  levels  of  the  laser  excitations  are  identical. 
Forinally,  the  ratio  of  tlie  signal  voltages  at  the  center  of  the  0^(1)  lines 
for  (1-0)  excitation  compared  to  (1-1)  excitation  is: 


^^s>i-i 

<^s^l-0  (X^Af^Q^)^^_Q^y[OH(v"  = 0)]^ 

The  ratio  of  signal  voltages  to  be  expected  for  (1-0)  and  (0-0)  excitation 

differs  because  the  upper  level  of  the  laser  induced  OH  transition  differs. 

This  ratio  must  account  then  for  electronic  quenching,  V-V  transfer  in  the 
2 + 24 

A H state,  rotational  redistribution  in  v'  = 1 and  v'  = 0,  and  the  wave- 
length dependent  filter  transmission  function.  For  our  purposes,  it  was 
easier  to  measure  the  integrated  effect  of  these  differences  and  apply  an 
empirically  determined  correction  factor,  Cp,  where 

'V, 


(b) 


C 


1-0 

'o-Oy 


(c) 
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Combining  expressions  (a),  (b)  and  (c)  leads  to  a numerical  value  for  Kj, 
for  use  in  the  following  expression 


m.  RESULTS  AND  ANALYSIS 


A typical  OH  excitation  fluorescence  spectrum  obtained  in  this 
experiment  is  shown  in  Fig.  6.  This  was  obtained  by  simultaneously  tun- 
ing the  laser  and  frequency  doubling  crystal  over  the  wavelength  interval 
shown  and  recording  the  output  of  the  Boxcar  Averager.  The  0^(1)  line 
was  used  for  quantitative  work  because  it  is  a strong  transition  with  a 
large  fractional  population  in  the  lower  level  and  it  is  well  isolated  from 
possibly  interfering  lines. 


Data  for  the  dependence  of  OH(v  = 0,  1)  density  on  the  flow  rate  of 
NO  are  shown  in  Fig.  7.  The  rate  constants  reported  in  this  paper  are 
derived  from  these  and  similar  data.  These  were  obtained  for  constant 
Inlet  N-atom  density  and  H2(v"  = 1)  density.  These  data  were  obtained  by 
changing  the  NO  gas  flow  rate  for  chronologically  sequential  data  points 
In  a random  way,  so  as  to  minimize  the  effect  of  systematic  trends. 


The  kinetics  of  the  O + H2  system  were  modeled  using  a computer 
27 

code  which  integrates  the  chemical  rate  equations  for  a constant  area 
streamtube  consistent  with  the  conservation  of  energy,  momentum,  and 
atomic  species.  Because  of  the  rapid  reactions  (3)  through  (9)  listed 
below,  the  densities  of  OH(v  = 0)  and  OH(v  = 1)  were  steady  state. 

OH(v  =0):  O + OH(0)— ► ©2  + H (3) 

H2  + OH(0)  H2O  + H (4) 

N + OH(0)  — *NO  + H (5) 
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I Fig.  6.  An  Excitation  Fluorescence  Spectrum  of  the  (0-0)  Band  in  OH. 

The  OH  was  generated  by  reaction  in  the  O + H2  system.  The 
lines  are  identified  according  to  the  nomenclature  of  Ref.  14. 

1 


I 

i 
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FLOW  RATE  RATIO,  ( 

Fig.  7.  Data  for  the  Density  of  OH(v  = 0,  1)  As  a Function  of  the 
Flow  Rate  of  NO  for  Constant  N-Atom  Flow  Rate 


OH(v=l);  O + OH(l)— + H (6a) 

O + OH(l)— >0  + OH(0)  (6b) 

H + OH(l)  — H + OH(0)  (7a) 

NO  + OH(l)— ►NO  + OH(0)  (8a) 

N + OH(l)  — ►NO  + H (9a) 

N + OH(l)— ► N + OH(0)  (9b) 


Figure  2 shows  why  some  of  the  reactions  (3)  through  (9)  become  unimpor- 
tant and  others  become  dominant  depending  on  the  magnitude  of  §.  The 
N-atom  reactions  are  negligible  for  §2  1.0;  the  O-atom  reactions  are 
unimportant  and  the  N-atom  reactions  become  dominant  for  5 « 1.  For 
OH(v  = 1)  the  quenching  by  NO  becomes  the  only  important  loss  process 
for  I i 8.  0.  For  these  reasons,  the  analysis  is  simplified  if  we  concen- 
trate on  restricted  ranges  of  values  of  § and  consider  other  §-value  regions 
separately. 

In  the  following,  steady  state  analyses  are  presented  first.  These 
lead  to  an  approximate  value  for  for  § » 1 and  give  physical  insight 
regarding  the  important  reactions  for  all  § values.  While  the  densities  of 
OH(v  = 0,  1)  are  steady  state,  the  densities  of  other  species  important  in 
the  analyses  are  not  steady  state.  Final  values  of  rate  constants,  there- 
fore, are  found  from  a computer  fit  to  the  data  for  [OH(v)]  versus  5* 
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where. 


k3[0]  + kjH2(0)] 

H + OH(2)  — ► H + OH(l)  (7b) 

NO  + OH(2)— ►NO  + OH(l)  (8b) 

H + NO^— ► OH(v)  + NO  (10) 

and  10a,  10b,  10c,  lOd  refer  to  the  formation  of 
OH(v  = 0,  1,  2,  3),  respectively,  and 
O + HNO  — ► OH(0)  + NO  (11a) 

O + HNO— ►OH(l)  + NO  (11b) 

The  NO^  in  reaction  (10)  is  formed  from  the  reaction 

O + NO  + M— ► NO^  + M 

The  HNO  in  reaction  (11)  is  formed  in 

H + NO  + M— ► HNO  + M 

Thus,  both  NO2  and  HNO  are  present  in  significant  quantities  only  for 
5»  1.  The  OH(2)  appearing  in  expression  (e)  comes  from  reaction  (10) 


Of  the  terms  in  the  numerator  of  expression  (e),  the  first  one  involv- 
ing ^2a,  largest  value.  The  second  term  involving  reaction  (10b) 

is  of  comparable  magnitude  only  for  § > 25.  The  third  and  fourth  terms 
are  smaller  still.  Thus,  to  first  approximation  we  can  ignore  the  second. 


third  and  fourth  terms  to  arrive  at  cin  approximate  expression  for  k,  . 


k_  « 
2a 


kgjNO][OH(v  = 1)] 


13 


[0][H2(v  = 1)] 

Glass*"'  has  recently  obtained  a value  for  k^,  . The  densities  of  H_(v  = 1) 

03.  C 

and  OH(v  = 1)  are  available  from  the  experiment  and  the  densities  of  O and 
NO  can  be  calculated  accurately  from  the  measured  NO  flow  rates  at  5 = 1 
and  at  the  condition  for  which  [OH(v  = 1)]  is  measured. 

A corresponding  simplification  of  expression  (f)  to  arrive  at  an 
expression  for  k^^^  is  not  possible.  This  occurs  because  for  a good  match 
between  experimental  data  and  computed  values  for  [OH(0)],  the  value  of 
k^^  must  be  made  so  small  as  to  be  negligible  in  comparison  with  the  sum 
of  the  other  terms  in  the  numerator  of  expression  (f).  Furthermore,  there 
is  no  region  in  the  range  of  §-values  (5  » 1)  where  one  of  the  remaining 
terms  in  the  numerator  of  (f)  is  clearly  dominant.  Thus,  a simple  expres- 
sion is  not  obtainable  from  which  any  one  of  the  rates  k_,  , k,.  , k^  , k_  , 

CD  xua  oa  (2i 

or  kg^  could  be  deduced  in  terms  of  the  reasonably  well  known  values^  for 


(g) 


kj  and  k^,  and  the  density  of  OH(0).  However,  since  values  are  available 

28,  29,  30 

in  the  literature  for  k^Q,  k^^ 

uppjer  limit  to  k2^  can  be  established. 


13,17.  13,  13  ,,  31 

^8a'  ^7a'  ^6b' 
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The  data  from  which  a value  of  was  deduced,  by  means  of 


expression  (g)  are  given  in  Table  II.  From  these  we  find  k_  “•  1.8  x 10 


14 


cc/molecule-sec  using  the  value  of  kg^  from  Ref.  13.  This  is  an  approx- 
imate value;  a more  accurate  value  is  obtained  from  the  computer  simula- 
tion which  accounts  for  the  processes  in  expression  (e)  which  were 
neglected  in  deriving  expression  (g). 


J , 


Table  II. 


r 

» 

Sfi 
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, Experimental  Data  Used  to  Find  (^2a^^8a^  from  Expression  (g) 


5 

[NO]^ 

[o] 

[H2(v  = 1)] 

[OH(v  = 1)] 

6.5 

2.08(14)^ 

2.41(13) 

1.93(14) 

2.96(9) 

1.32(-4) 

6.  5 

2.  08(14) 

2.41(13) 

1.  93(14) 

2.65(9) 

1. 19(-4) 

9.3 

3.  13(14) 

2. 11(13) 

1.93(14) 

3.34(9) 

2.  57(-4) 

9.3 

3.  13(14) 

2. 11(13) 

1.93(14) 

2.88(9) 

2.21(-4) 

11.8 

7.  05(14) 

2.74(13) 

1.  93(14) 

2.55(9) 

3.40(.4) 

11.8 

7.05(14) 

2.74(13) 

1.93(14) 

1.57(9) 

2.09(-4) 

13.5 

4.75(14) 

1.87(13) 

1.  93(14) 

2.43(9) 

3.20(-4) 

13.5 

4. 75(14) 

1.87(13) 

1.93(14) 

2.12(9) 

2.79(-4) 

17.6 

1.08(15) 

1.93(13) 

1.99(14) 

1.18(9) 

3.32(-4) 

17.6 

1.08(15) 

1.93(13) 

1.99(14) 

1.18(9) 

3.  32(-4) 

24.0 

8.73(14) 

1.39(13) 

1.93(14) 

6.90(8) 

2.25(-4) 

24.0 

8.73(14) 

1.39(13) 

1.93(14) 

6.37(8) 

2.07(-4) 

Averag 

e value: 

<^2a'^^8a^  " 

48  ± .76)  X 10“* 

*Densitie8  are  expressed  in  molecules/cc. 

^By  definition:  2.08(14)  k 2.  08  x 10^^  in  this  table. 
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b.  § » 1 


The  steady  state  expressions  In  this  region  are; 

[OH(l)]  = 

k2^[0][H  (1)] 

2 

" ‘'6b>  ^ " *'l4b> 

and 

[OH(0)]  = 

[O]  (kitH^(O))  4 k^jH^d)])  4 [OH(l))(k^JOl  4 k^jH]  4 

k3[0]  + yH^CO)] 

where, 

H2(0)  + 0H(1)->H20  + H (14a) 

H2(0)  + OH(l)  — H^(l)  + OH(0)  (14b) 

Expression  (h)  shows  that  there  is  only  one  source  of  OH(v  = 1), 

1.  e.  , reaction  (2a),  but  now  there  are  several  Important  loss  processes, 

no  one  of  which  Is  clearly  dominant.  The  reactions  (14a)  and  (14b)  for 

32 

which  Spencer  et  al.  have  determined  upper  limits  only,  are  not  as 
Important  as  reactions  (6)  and  (7a).  The  H-atoms  are  produced  mostly  by 
reaction  (2)  with  minor  contributions  from  (3)  and  (4),  and  we  have  no 
experimental  measurement  of  [H],  This  makes  the  determination  of  a rate 
constant  a nonlinear  problem  which  is  best  handled  by  the  computer  simu- 
lation. 

With  regard  to  OH(0),  the  situation  is  qualitatively  as  stated  pre- 
viously for  § » 1.  In  order  to  match  the  data,  the  value  of  k,,,  must  be 

2b 
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so  small  as  to  make  reaction  (2b)  a negligible  source  of  OH(0).  The  major 
sources  of  OH(0)  are  now  the  physical  deactivation  of  OH(l)  by  H and  O. 

It  is  true,  though,  that  the  total  deactivation  rate  of  OH(v  = 1)  is  smaller 
for  § « 1 than  it  is  for  § » 1.  This  is  because  the  NO  density  is  greatly 
reduced  (see  Fig.  2)  while  the  densities  of  the  other  major  deactivating 
species  stay  approxiiTiately  the  same.  This  means  that  an  experimentally 
derived  upper  limit  to  is  bet.t  determined  from  data  where  § « 1.  since 
for  all  other  conditions  reaction  2b  is  a correspondingly  smaller  fraction 
of  the  total  rate  of  formation  of  OH(0). 


c.  I « 1 

In  this  region  it  is  found  that  the  unknown  rate  constants  for  proces- 
ses (9a)  and  (9b)  must  play  an  important  role.  Therefore,  it  is  not  possi- 
ble to  learn  anything  further  regarding  the  reactions  of  interest  by  analy- 
zing data  in  this  region.  However,  using  the  values  for  k^^  and  k^j^ 
determined  for  C ^ 1 U is  possible  to  deduce  values  for  k^^  and  k^^^  from 
the  data  for  % < 1.  This  is  to  be  the  subject  of  another  paper  to  be  pub- 

lished  in  the  near  future. 


d.  Final  Rate  Constants 

Extensive  computer  calculations  were  carried  out  to  determine  the 
best  values  lot  and  the  upper  Unait  to  k^^.  The  value  ol  k^^,  which 
brings  calculated  and  measured  values  ol  [OH(l)l  into  agreement,  was 
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determLned  for  ajl  data  for  which  5^1.  This  includes  more  data  than  are 
shown  in  Table  II.  The  average  of  these  values  was  k,  = (1. 0 ± 0.  4)  x 10” 
cc/molecule- sec.  The  uncertainty  is  an  experimental  uncertainty  and 
represents  one  standard  deviation.  The  computed  densities  of  OH(v  = 0,  1) 
for  this  value  of  k are  shown  in  Fig.  7 as  the  solid  lines.  For  this 

di 

calculation  the  value  of  k-,  was  set  to  zero.  The  values  used  for  all 

CD 

other  important  rates  in  this  computation  are  listed  in  Appendix  B. 


-14 

With  k^^  = 1.  0 X 10  , computations  were  done  for  a range  of 

values  of  k^,  . These  were  done  with  the  values  of  rate  constants  for 
cD 

other  reactions  which  affect  [OH{0)]  chosen  to  minimize  their  contribution 
to  [OH(0)],  consistent  with  forming  an  upper  limit  to  k^^.  The  value  of 
^2b  brought  about  agreement  between  calculated  and  measured 

[OH(0)]  was  determined  for  each  data  point  in  the  range  0.  9 ^ ^ 2.  0. 

The  average  value  was  k^^  ^ 4.  7 x 10  cc/molecule-sec.  This  must  be 
taken  as  an  upper  limit  to  the  true  value  since  Fig.  7 shows  that  computa- 
tions with  k2|j  = 0.  0 give  a reasonably  good  fit  to  the  data. 

e.  Sensitivity  Analysis 

The  values  of  the  rate  constants  reported  here  are,  of  course,  sen- 
sitive to  uncertainties  in  the  rates  of  certain  other  reactions  as  used  in  the 
computer  simulations. 


14 
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" OHIv=0) 

EXPERIMENTAL  DATA 
o OH  (v=OI 
X OH  (v=ll 

- COMPUTED  DENSITIES 


Fig.  8. 


FLOW  RATE  RATIO, ( 


Comparison  of  the  Computer -Calculated  Densities  of  OH(v  = 0,  1) 
With  the  Experimental  Data  Shown  in  Fig.  7.  For  this  calculation, 
k,  =1.0x10"^'^  cc/molecule -sec  and  k,,  = 0.  0. 


The  most  important  sources  of  OH(v)  in  these  experiments  other 
than  reactions  (2)  are  the  following  reaction  sequences: 

O + NO  + M — ►NO^  + M U2) 

followed  by 

H + NO^ — ►OH(v  s 3)  + NO  (10) 

and 

H + NO  + M — ► HNO  + M (13) 

followed  by 

O + HNO — ►OH(v  s 5)  + NO  (H) 

Experimentally  measured  values  of  reaction  rate  have  been  reported 

for  all  of  these  except  reaction  (11).  There  has  been  extensive  recent 

work  on  each  of  the  other  three  reactions.  These  three  were  also  reviewed 

by  Baulch  et  al.  Baulch's  suggested  uncertainty  in  reaction  (12)  is 

± 20%  and  measurements  have  been  reported  for  M = N^  and  M = the 

principal  constituents  in  the  experiments  reported  here.  The  overall  reac- 
* ^ 28  29  30 

tion  rate  for  reaction  (10)  has  been  the  subject  of  three  recent  studies. 

These  three  yield  values  for  which  are  in  substantial  agreement  with  one 

another  but  which  differ  by  a factor  of  about  2.  5 from  the  previously  held 

value.  It  seems  certain  that  the  more  recent  values  are  to  be  preferred. 

When  the  data  from  Refs.  28  29.  and  30  are  weighted  equally  the  value, 

k = (1.22  ± .08)  X 10‘^°  is  obtained^®  for  T = 298 “k.  It  is  also  true  that 

the  vibrational  branching  ratios  for  reaction  (10).  although  they  have  recently 

received  much  attention/ ^ are  still  not  accurately  known.  Spencer 

and  Glass^^  have  used  Bemand's  value^^  k^^  = 1.  I X lO’^®  in  interpreting 

their  measurements  of  the  quenching  rates  for 
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NO  + OH(v)— ► NO  + OH(v  - 1)  (8) 

H + OH(v)— ►H  + OH(v  - 1)  (7) 

Spencer  and  Glass  have  presented  their  data  for  these  rates  In  such  a way 

that  new  values  can  be  derived  for  different  values  of  the  overall  rate  and 

branching  ratios  of  reaction  (10).  We  have  reinterpreted  Spencer  and 

Glass's  results  for  the  rates  of  (7)  and  (8)  for  kjQ  = 4.  8 X 10 

cc/molecule-sec  and  for  a set  of  vibrational  branching  ratios  consistent 

17  X 6 

with  the  work  of  Silver  et  al.  and  Polanyi  and  Sloan.  The  relative 

rates  are  1. 0/1.  3/1. 3/0.  13  for  V = 0 to  3,  respectively.  In  assessing  the 

overall  uncertainty  in  k^^.  the  values  k^^  = 1.  3 x 10  the  ratios 

k,«./k,  » /k,_,  = 1.0/1.3/1.3/0.13,  respectively  and  Spencer's  reported 
lUb  lUc  lUa 

values  for  k.,  and  kj,  were  used  in  determining  a lower  limit,  while 

kio=  4.8  X 10-",  ki0a'''l0b^‘'l0c^‘‘l0d  ' 1- 0/0- ''/O.  3/0.  03  and  appro- 
priately  modified  values  of  k.j,  and  kg  were  used  in  determining  an  upper 
limit. 

Baulch  et  al.  suggest  an  uncertainty  of  ± 50%  for  the  rate  of  (13) 
with  M = No  values  seem  to  have  been  reported  for  M = N^.  We  shall 

assume  that  the  rate  for  M = N^  is  the  same  as  for  M = for  purposes  of 
this  uncertainty  analysis. 

35  “11 

While  a lower  limit  has  been  placed^  on  k^  ^ ( 2 1.8  x 10  ), 

nothing  is  known  about  the  branching  ratios  for  forming  OH(v).  With 
kjj  = 1.8  X 10'^^,  then  reaction  (13)  is  rate  limiting  and  larger  values 
of  kj  j will  have  no  effect  on  the  calculated  OH  density. 
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For  our  purposes  then,  it  is  quite  sufficient  to  assume  that  = 0.0 

and  k, = 1.  8 x lO”^^  cc/molecule-sec.  , i.  e,  , that  only  OH(v  = 1)  is  formed 
lib 

when  a lower  limit  to  is  being  evaluated;  conversely  we  assume  that 

k,  , = k, = 0 when  an  upper  limit  is  being  evaluated. 

11a  11b 

One  other  issue  needs  to  be  addressed  with  regard  to  these  possibly 
interfering  reactions  (10)  through  (13).  H-atoms  are  formed  by  thermal 
dissociation  of  on  the  hot  tungsten  filament.  They  also  are  formed  by 
the  reactions  of  interest, 

O + ►OH  + H. 

If  one  calculates  the  equilibrium  dissociation  level  at  the  measured  tem- 
perature of  the  tungsten  filament  one  finds  that  the  H-atom  flow  rate  con- 
tributed by  the  filament  is  roughly  comparable  to  that  formed  by  subsequent 
reaction  after  the  reactants  are  mixed.  This  means  that  uncertainties  in 
the  density  of  H-atoms  assumed  to  be  coming  from  the  filament  might  be 
important  in  determining  the  best  value  for  k2j^,  since  [H]  enters  the  pic- 
ture through  both  reactions  (10)  and  (13)  as  well  as  reaction  (7).  Attempts 
were  made  early  in  this  project  to  infer  the  H-atom  flow  rate  from  the 
filament  in  separate  experiments.  These  consisted  of  allowing  the  hydro- 
gen coming  from  the  hot  filament  to  mix  with  a stream  of  carrying  a 
known  amount  of  NO2  gas.  The  products  of  reaction  were  examined  for 
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> OH  molecules  by  laser  induced  fluorescence.  The  absolute  density  of 

I *• 

1 pir‘ 

[ OH(0)  was  used  in  conjunction  with  the  known  NO^  flow  rate  to  infer  an 

I V,i  initial  H-atom  density.  The  result  of  these  measurements  was  that  the 

inferred  value  of  [H]  was  at,  or  slightly  below  the  equilibrium  value  over 

' * O O 

the  filament  temperature  range  from  800  C to  1150  C, 


Accordingly,  the  simulation  calculations  were  performed  with  an 
H-atom  flow  rate  determined  by  equilibrium  dissociation  at  the  filament 
temperature. 


Computer  calculations  were  carried  out  in  which  the  rates  of  all 
important  reactions  were  allowed  to  vary  over  the  ranges  indicated  in 
Appendix  B.  The  process  of  determining  a value  for  for  each  data 
point  was  repeated  to  find  upper  and  lower  limits  to  k^^  consistent  with 
these  rate  constant  uncertainties.  The  results  of  this  analysis  are; 

k.,  = (1.0  ^ n"/ ) X lO”^^  cc/molecule-sec.  It  is  seen  that  the  experi- 

* U • o 

- 14 

mental  uncertainty  of  (±  0.4)  X 10  is  not  greatly  different  from  the 
uncertainty  arising  from  uncertain  rate  constants.  Thus,  our  best  esti- 


mate of  ko  is; 
2a 


'^Za  “ ^ 0*6^  ^ 10”^^  cc/molecule-sec, 


With  regard  to  the  upper  limit  to  the  values  for  the  rates  of 

competing  reactions  were  already  set  to  their  respective  uncertainty 
limits  for  the  computations  described  in  Section  d.  Hence,  the  reported 
upper  limit  to  ^21)  ^^kes  into  consideration  these  other  uncertainties. 
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; 


The  reported  value  for  and  upper  limit  to  can  be  summed 

3 

to  obtain  the  total  reaction  rate  of  0(  P)  with  H2(v  = 1).  Thus, 

“ ^2a  ^ ^2b  cc /molecule- sec. 


JT'’! 


The  question  naturally  arises  whether  there  might  be  contributions 
to  the  observed  reaction  rate  from  vibrationad  levels  in  H2  higher  than 
V = 1.  Vibrationally  excited  is  generated  by  a thermal  source  in  this 
experiment  and  there  is  every  reason  to  expect  that  higher  vibrational 
states  are  represented  in  the  flow.  The  vibrational  temperature 
implied  by  measured  values  of  [H2(v  = 1)]  is  T^  » 1100  K.  Even  if  the 
original  source  were  not  thermal  there  would  be  a rapid  redistribution  of 
vibrational  energy  by  nearly  resonant  V-V  exchange  between  H2  molecules. 

While  it  cannot  be  rigorously  proven  that  higher  vibrational  levels 
(v  St  2)  are  not  contributing  to  the  observed  OH  production,  the  following 
strongly  inferential  observation  has  been  made.  Data  were  obtained  for 
the  dependence  of  OH(v  si)  on  the  density  of  H2(v  = 1).  These  data  will 
be  the  subject  of  a future  publication.  If  H2(v  2 2)  were  significant  the 
dependence  of  [OH(v)]  on  [H2(1)]  would  be  expected  to  be  of  order  larger 
than  one.  In  fact,  the  data  show  that  [OH(v)]  is  of  order  one  or  less  in 
[H2(1)].  Thus,  it  is  felt  that  H2(v  a 2)  does  not  contribute  to  the  observed 
OH(v)  and  that  the  measured  rate  constant  refers  strictly  to  H2(v  = 1) 
reactant. 

The  largest  possible  value  of  k2  consistent  with  the  uncertainties 
— - 14  ' 

determined  above  is  (k-)  = 1.9  X 10  . This  is  smaller  by  a comfor- 


The  reported  value  represents  an  increase  owing  to  vibrational  excita- 
tion  in  hydrogen  by  a factor  of  » 2.  6 X 10^  over  the  rate  of  reaction  (1), 

O 

although  the  rate  of  reaction  (1)  has  never  been  measured  at  T 300  K. 

The  acceleration  of  a reaction  process  by  vibrational  energy  in  one  of  the 

36 

reactants  is  characterized  by  some  authors  by  the  fraction,  a,  of  the  vibra- 
tional quantum,  which  is  effective  in  reducing  the  activation  energy 

for  that  reaction.  From  this  approach  the  ratio  k^/k  of  vibratlonally 

accelerated  rate  to  ground  state  rate  constants  is 

* 

k /k  = exp(o  Eyjg/RT). 

For  the  value  of  k^  reported  herein,  it  is  found  that  or  « 0.  43. 


There  are  two  classical  trajectory  treatments  of  reaction  (2)  in  the 

g 

literature  with  different  approaches.  Johnson  and  Winter  did  a collinear 
OHH  analysis  using  an  LEPS  potential  surface.  The  barrier  height  was 


first  adjusted  to  bring  about  agreement  between  the  calculated  and  meas- 

5 4 9 

ured  ' rate  constants  for  the  fully  ground  state  reaction.  Whitlock  et  al. 

generated  a diatomics-in-molecules  potential  surface  which  confirmed 

3 

that  collinear  hydrogen  abstraction  was  the  dominant  0(  P)  + reaction 

path.  Johnson  and  Winter  report  an  expression  for  the  rate  constant  for 

- 14 

reaction  (2)  of  k^  = 4.  65T  x 10  exp(-1868/T).  At  room  temperature 

-14 

this  gives  k2  = 2.  9 X 10  , a value  in  remarkably  good  agreement  with 

the  experimental  value  reported  herein.  Whitlock  et  al.  report  a range  of 

values  for  k2  of  1.  9 X 10  s k2  ^ 8.  3 x 10  at  T = 300  K.  Johnson  and 

Winter  also  report  an  expression  for  the  branching  ratio,  F s 1^2a^^2b 

- 14 

r = 2.  3 exp(196/T).  This  leads  to  a value  o[  = 2.27  x 10  , again  in 


reasonably  good  agreement  with  experiment.  The  theoretically  deter- 

- 15 

mined  value  for  i.  e. , k2j^  = 5.  15  X 10  at  room  temperature  is 

- 1 5 

very  close  to  the  experimental  upper  limit  of  s 4.  7 x 10 
cc /molecule- sec. 


The  good  agreement  between  the  result  of  Johnson  and  Winter  and 
the  present  experimental  result  does  not  necessarily  attest  to  the  accuracy 
of  their  potential  surface  because  of  the  uncertainties  introduced  by  a clas- 
sical treatment  of  vibrations  and  because  tunneling  is  not  modeled.  The 
results  of  Ref.  9 confirm  the  predominance  of  the  collinear  abstraction 
mode  but  give  a considerably  different  predicted  rate  constant. 

The  present  understanding  of  the  effect  of  vibrational  energy  on  the 
reaction  of  O-atoms  with  is  summarized  graphically  in  Fig.  9.  It  is 
seen  that  the  theoretical  expression  for  k^  from  Ref.  8 lies  slightly  above 
the  uncertainty  range  of  the  experimental  measurement  reported  herein. 

It  is  of  interest  to  examine  the  implications  of  the  rate  of  reaction 
(2)  for  the  measurement  of  the  rate  of  reaction  (1)  at  high  temperatures. 

At  sufficiently  high  temperatures,  assuming  thermal  equilibrium  between 
translational  and  vibrational  modes,  there  will  be  enough  H2(v  = 1)  present 
to  affect  the  apparent  rate  of  reaction  (1).  Indeed,  recent  measurements^ 
of  kj  at  high  temperatures(1400“K  s T s 1900“K)  give  values  which  are 
significantly  larger  than  predicted  by  a simple  Arrhenius  expression 
which  fits  all  of  the  d ita  in  the  temperature  range  900  a T 4 400*K.  The 
possibility  that  this  upward  curvature  at  high  temperatures  might  be  the 
result  of  reaction  (2)  was  suggested  by  Gardiner.  It  should  be  pointed 
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out  here  that  non-Arrhenius  behavior  is  not  a particularly  surprising 
result  even  for  a fully  ground  state  reaction.  It  is  the  purpose  of  this 
exercise  to  establish  only  whether  or  not  the  vibrationally  accelerated 
reaction  might  be  responsible  for  the  observed  curvature. 

We  adopt  the  Arrhenius  expression  given  by  Johnson  and  Winter  for 

O 

k^,  which  is  in  good  agreement  with  our  experimental  result  at  T = 302  K, 
and  calculate  the  total  rate  of  reaction  k , including  the  effect  of  H (v  = 1) 

t M 

according  to  the  following  expression: 

kj[H2(v  = 0)]  + k2[H2(v  = 1)] 

‘ i:[H,(v)] 


For  kj  we  adopt  the  consensus  expression  presented  by  Baulch  et  al. 

2 

The  result  is  shown  in  Fig.  10  along  with  the  data  of  Schott  et  al.  , the 
results  of  a transition  state  theory(t.  s.t.  )-BEBO  calculation  also  presented 
by  Schott,  and  the  Baulch  expression  for  k^.  It  is  seen  that  k^,  according 
to  expression  (j)  above,  does  not  match  the  data  of  Schott  et  al.  as  well  as 
do  the  t.  s.  t.  calculations.  The  differences  are  small,  however,  and  a final 
resolution  must  await  data  for  the  value  of  k2  in  the  temperature  range 
T ~ 2000  "k. 


The  observation,  reported  herein  for  the  first  time,  of  the  existence 
of  vibrational  excitation  in  H2  when  is  mixed  with  active  nitrogen  may 
have  implications  for  previous  measurements  for  T < 400 °K  of  the  fully 
ground  state  reaction  rate  constant,  i.  e.  , reaction  (1).  There  have  been 
at  least  three  reports  ’ ’ of  the  value  of  this  rate  constant  for  which  the 
unknown  existence  of  vibrationally  excited  hydrogen  produced  by  energy 


L 


CONSTANT,  k,  (cc/molecule-seci 


exchange  from  active  could  have  been  a factor.  Among  these  three 

reports  there  are  significant  differences  in  the  reported  rate  constant 

values  at  low  temperatures  where  the  effect  of  undetected  vibrationally 

excited  would  be  most  pronounced.  We  have  obtained  measurements 

for  the  rate  of  reaction  (1)  at  a lower  temperature  than  previously  reported 

by  any  investigators.  These  measurements  were  obtained  with  careful 

monitoring  and  accounting  for  the  effects  of  possible  H2(v  = 1).  These 

38 

results  will  be  reported  separately. 
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Appendix  A 


Reactions  Used  for  the  Calculation  of  Species  Concentrations 

in  the  N + NO  System 


Included  here  are  the  specifications  and  rate  coefficients  for  all  of 
the  chemical  reactions  used  to  calculate  species  concentrations  resulting 
from  the  titration  of  active  nitrogen  with  NO.  The  rate  coefficients  are 
represented  by  three  coefficients  a,  b,  c such  that  the  rate  constant  k is 
given  by 

k = aT^exp(c/T),  cc/mole-sec 


where  T is  in  degrees  Kelvin.  The 
from  that  used  in  the  text. 

Reaction 

(1)  N + NO— ► N^  + O 

(2)  20  + M— ► ©2  + M 

(3)  0 + 0^  + M— ►Og  + M 

(4)  O + NO  + M — ► NO2  + M 

(5)  O + N + M — ► NO  + M 

(6)  2N  + M— ► N2  + M 

(7)  O + O3— ► 2O2 

(8)  N + O2— ► NO  + O 

(9)  O + NO2  — ► O2  + NO 

(10)  O + NO— ► N02(+  hv) 

(11)  NO  + O^— ► NO2  + ©2 

(12)  N+NO2— ►N2O  + O 


numbering  of  the  reactions  is  different 


a 

b 

c 

3.  1 

X 

10'' 

0 

-168, 

1.4 

X 

lo'S 

-1.0 

-171 

1.7 

X 

10^3 

0 

+ 1057, 

1.05 

X 

10^^ 

0 

+941, 

6.  4 

X 

lo'^ 

-0.  5 

0 

2.8 

X 

10^4 

0 

+503 

1.2 

X 

10^^ 

0 

-2416 

6.  4 

X 

10^ 

+ 1.  0 

-3145 

1.0 

X 

10^3 

0 

-302 

3.4 

X 

10^2 

-2,0 

0 

5.7 

X 

lo'i 

0 

-1309 

9.03 

X 

10'' 

0 

0 

-59. 


h 


In  addition,  the  removal  of  O at  the  wall  was  simulated  by  specifying  a 
first  order  rate  constant,  = -3.4  sec  , for  use  in  the  expression; 


wall 


k [O] 
w 
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Appendix  B 

Reactions  Used  for  Computer  Calculations  in  the  O + System 

Assembled  here  are  the  ZO  most  important  chemical  reactions  and 

rate  coefficients  used  in  the  computer  analyses  described  in  the  text  for 

the  O + system,  A total  of  102  reactions  were  used  in  each  of  the 

computations  but  only  20  of  these  proved  to  be  important.  Only  those 

reactions  which  significantly  affected  the  concentration  of  any  species  are 

included  in  the  following  list.  In  particular,  many  reactions  involving  the 

production  and  destruction  of  HO^,  HNO^,  HNO^,  and  N^O  were 

39 

used  in  the  computation  but  do  not  appear  below  because  none  of  them 
were  significant  for  calculated  concentrations  of  OH(0,  1).  The  nomen- 
clature is  the  same  as  used  in  Appendix  A but  the  numbering  of  the  reac- 
tions is  the  same  as  that  used  in  the  text. 


Reaction 

a 

b 

c 

Uncertainty 

(1) 

O + H2(0)— *OH(0)  + H 

1.8 

X 

lo'O 

+1.0 

-4480. 

- 

(3) 

O + OH(0)— ^02  + H 

2.3 

X 

10^^ 

0 

0 

±30% 

(4) 

H2(0)  + OH(0)— ► H2O  + H 

2.2 

X 

10^^ 

0 

-2590. 

±20% 

(5) 

N + OH(0)  — ► NO  + H 

3.2 

X 

10^3 

0 

0 

- 

(6a) 

0 + OH(l)— *>02  + H 

6.32 

X 

10^^ 

0 

0 

±50% 

(6b) 

0 + OH(l)— ► 0 + OH(0) 

8.73 

X 

0 

0 

±20% 

(7a) 

H + OH(l)— ►H  + OH(0) 

1.63 

X 

10'^ 

0 

0 

+0% 

-56% 

(7b) 

H + OH(2)— t-H  + 0H(1) 

1.99  X 

10^^ 

0 

0 

+ 9% 
-56% 
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Reaction 

a 

b 

c 

Uncertainty 

(8a) 

NO  + OH(l)— ►NO  + OH(0) 

4.  52 

X 

1013 

0 

0 

+0% 

-56% 

(8b) 

NO  + OH(2)— ► NO  + OH(l) 

7.23 

X 

1013 

0 

0 

+ 9% 

-56% 

(10a) 

H + NO2— ►OH(O)  + NO 

3.71 

X 

10^3 

0 

0 

+6% 

-80% 

(10b) 

H + NO^— ♦OH(l)  + NO 

2.  59 

X 

10^^ 

0 

0 

+ 6% 

-60% 

(10c) 

H + NO2— ►OH(2)  + NO 

1.04 

X 

10^^ 

0 

0 

+6% 

-60% 

(11a) 

0 + HNO— ►NO  + OH(0) 

6.  0 

X 

10^2 

0 

0 

+0% 

-100% 

(11b) 

0 + HNO— ►NO  + OH(l) 

6.0 

X 

10^2 

0 

0 

+ 0% 
-100% 

(12a) 

0 + NO  + N^— ► NO2  + N^ 

1.  54 

X 

10^^ 

0 

+940. 

±20% 

(12b) 

0 + NO  + — ► NO^  + 

1.82 

X 

10^^ 

0 

+940. 

±20% 

(13) 

H + NO  + M — ►HNO  + M 

5.4 

X 

10^5 

0 

+300. 

±50% 

(14a) 

H2(0)  + OH(l)— H^O  + H 

6.0 

X 

10^ 

0 

0 

+0% 

-30% 

(14b) 

H2(0)  + OH(l)— H2(1)  + OH(0) 

6.  0 

X 

10^ 

0 

0 

+0% 

-100% 

In  addition,  removal  at  the  wall  was  handled  by  assigning  the  follow- 
ing first  order  rate  constants: 


Species 

k , sec 
w 

0 

-3.4 

H 

-4.0 

H^d) 

-46.0 

OH(0) 

-30.0 

OH(l) 

-80.0 
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THE  IVAN  A.  GETTING  LABORATORIES 

The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 
fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 


matrix  composites  and  new  forms  ^ carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue-induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory;  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explo8ions> 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems. 

THE  AEROSPACE  CORPORATION 
El  Segundo,  California 
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